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Abstract 

A /3-beam with maximum 7 = 150 (for ®He ions) or 7 = 250 (for ^®Ne) could be achieved 
at the CERN-SPS. We study the sensitivity to ^13 and 6 of such a beam as function of 7 , 
optimizing with the baseline constrained to CERN-Prejus (130 km), and also with simultane¬ 
ous variation of the baseline. These results are compared to the standard scenario previously 
considered, with lower 7 = 60/100, and also with a higher 7 ~ 350 option that requires a 
more powerful accelerator. Although higher 7 is better, loss of sensitivity to @13 and 5 is most 
pronounced for 7 below 100 . 
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1. Introduction 


Results from atmospheric solar reactor j^I and long-baseline neutrino experiments 
in recent years can be economically accommodated in the Standard Model (SM) with neutrino 
masses and a three-neutrino mixing matrix [Sj. In this case, the lepton sector of the SM 
closely resembles that of the quarks and there are new physical parameters measurable at low 
energies: the three neutrino masses, rrii {i = 1, 2, 3), three mixing angles, 9ij, {i ^ j = 1, 2, 3), 
and a CP-violating phase, 5. In contrast with the quark sector, two additional phases could 
be present if neutrinos are Majorana. Of these new parameters, present experiments have 
determined just two neutrino mass-square differences and two mixing angles: (|Am 23 | — 
2.2 X 10-3 eV2, 023 ~ 45°) which mostly drive the atmospheric oscillation and {Ami 2 — 
8 X 10-5 eV^, 012 ~ 32°) which mostly drive the solar one. The third angle, 0i3, as well 
as the CP-violating phases (5, and possible Majorana phases) remain undetermined. Only 
an upper limit 0i3 < 12° is known. Another essential piece of information needed to clarify 
the low-energy structure of the lepton flavor sector of the SM is the neutrino mass hierarchy 
and the absolute neutrino mass scale. The former is related to the sign of the largest mass- 
square difference which determines if the spectrum is hierarchical (if the two most 

degenerate neutrinos are lighter than the third one) or degenerate (if they are heavier). 

Measurement of some of these parameters may be possible in high-precision neutrino- 
oscillation experiments. A number of experimental setups to significantly improve on present 
sensitivity to 0i3, 5 and the sign of Am^^ have been discussed in the literature: neutrino 
factories (neutrino beams from boosted-muon decays) [HI [Z1 El; superbeams (very intense 
conventional neutrino beams) 0 IIOL IllL I12| improved reactor experiments m and more 
recently /3-beams (neutrinos from boosted-ion decays) [m [isi. These are quite different 
in terms of systematics but all face a fundamental problem which limits the reach of each 
individual experiment significantly, namely the correlations and degeneracies between param¬ 
eters [ini-PH; 6*13 and 6 must be measured simultaneously, and other oscillation parameters 
are not known with perfect precision. 

To resolve these degeneracies it is important to measure as many independent channels 
as possible and to exploit the energy and/or baseline dependence of the oscillation signals 
and matter effects in neutrino propagation. In many cases, the best way to do this is by 
combining different experiments; indeed the synergies between some combinations of the 
setups mentioned above have been shown to be considerable. 

The neutrino factory provides ultimate sensitivity to leptonic CP violation, and thus rep¬ 
resents the last step on a long-term road map to reveal the lepton-flavor sector of the SM. 
Recently it was shown that a /3-beam running at a higher 7 than previously considered (and 
longer baselines), in combination with a massive water detector, can reach sensitivity to lep¬ 
tonic CP-violation and sign(Am 23 ) that competes with a neutrino factory’s. The optimal 
setup among those considered in uni was a /3-beam with 7 = 350/580 for ^He and 
isotopes respectively and a baseline L ~ 730 km. If constructed at CERN, this beam would 
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require a refurbished SPS or an acceleration scheme utilizing the LHC - implying snbstantial 
R&D effort in either case. 

This paper considers instead the possibility of using the existing CERN-SPS np to its 
maximum power, allowing a beam with 7 = 150 (250) for ®He (^®Ne) ions (some preliminary 
results of this study were presented in m)- The design of this /3-beam is essentially as 
described in m 

The advantages of increasing the 7 factor discnssed in [ 2 n] also apply in this case. The 
oscillation signals grow at least linearly with the 7 factor, therefore the highest 7 possible is 
preferred in principle, if the baseline is adjusted appropriately. Furthermore when the energy 
is well above the Fermi momentum of the target nuclei, energy dependence of the oscillation 
signals is very effective in resolving parameter degeneracies. In practice there are two caveats 
to this rule. First, water Cherenkov detectors are best suited for quasi-elastic (QE) reactions, 
where the neutrino energy can be kinematically reconstructed. Therefore sensitivity improves 
with 7 only until the inelastic cross section begins to dominate; we will show that this occurs 
for 7 > 400. The second concern is background, since NC single-pion production can mimic 
the the appearance signal; it is demonstrated in m and confirmed here that this background 
is manageable, even for 7 > 100 . 

More concretely the purpose of this paper is two-fold. First, optimization of a CERN-SPS 
/3-beam by answering the following questions: 

• Assuming an underground laboratory at Frejus with a megaton water Cherenkov de¬ 
tector, what is the optimal 7 with the existing CERN-SPS? 

• For the maximum 7 achievable with the CERN-SPS, what is the optimal /3-beam base¬ 
line? 

• Is there any physics advantage to varying the 7 ratio for ®He and ^®Ne, i.e. a ratio differ¬ 
ent from 7 i 8 Ne/ 7 ®He = 1-67 (which allows both beams to circulate simultaneously)? [ 2 SI 

Second, comparing the performance of the following set-ups: 

• Setup I: L = 130 km (CERN-Frejus) at the optimal 7 accessible to the CERN-SPS. 

• Setup II: 7 = 150 at the optimal baseline. 

• Setup III: 7 = 350 at L = 730 km, which is a symmetric version of the configuration 
considered in |25j . To accelerate the ions would require either a refurbished SPS (with 
superconducting magnets) or a more powerful accelerator, such as the Tevatron or LHC. 

In all cases an intensity of 2.9 x 10^® ®He and 1.1 x 10^^ ^®Ne decays per year |27| and an 
integrated luminosity corresponding to 10 years are assumed. Although these luminosities 
have been estimated for simultaneous ion circulation (fixing the ratio of 7 ’s to 1.67) refer¬ 
ence m argues they are achievable even if the ions circulate separately at the same 7 , by 
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Figure 1: (solid) and Ue (dashed) fluxes as a function of the neutrino energy at L = 300 km 

for the maximum acceleration of the ®He ('y = 150j and (7 = 250) at the CERN-SPS. 


injecting more bunches. While these intensities are realistic for the CERN-SPS, the same has 
not been demonstrated for other accelerators like the Tevatron or LHC. The far detector is 
a Super-Kamiokande-like water Cherenkov design, with fiducial mass 440 kton. 

The paper is organized as follows. Section 2 shows expected fluxes and event rates for the 
maximum 7 achievable at the CERN-SPS. Section 3 describes the performance of a large 
water Cherenkov detector for the appearance and disappearance signals and estimates the 
atmospheric background, an important constraint in design of the bunch length. Section 4 
deals with optimizations needed to define setups I and II and Section 5 compares the physics 
reach of the three emergent reference setups. Section 6 discusses our outlook and conclusions. 


2. Neutrino fluxes and rates 

Pigure^shows the fluxes for the maximum acceleration of the ions at the CERN-SPS: 7 = 150 
for ®He and 7 = 250 for ^®Ne at L = 300 km. Table 1 shows the rate of charged-current 
interactions expected per kiloton in one year. 
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7 

L{km) 

Pe CC 

I/e CC 

{E,){GeV) 

150/250 

300 

22.8 

115.6 

0.58/0.94 


Table 1: Number of charged-current events per kton-year, in the absence of oscillation, for 
the maximum acceleration of ®He and at the CERN-SPS. The average neutrino energy 
is also shown. 


3. Measurements at a /?-beam 

The parameters @13 and 8 are best studied by probing the appearance channels for neutrino 
oscillation in the atmospheric energy range: golden Vg) and silver (ur Vg) m 

channels have been identihed. In the setups considered here, neutrino energies are below r 
threshold, therefore only the golden channel is available. 

The disappearance transition Vg Vg can also be measured. This is an important comple¬ 
ment to the golden channel measurement, because the intrinsic degeneracy HZl in the golden 
measurement can be resolved: the disappearance measurement depends on @ 13 , but not on 
8. The synergy between the appearance and disappearance channels for a /3-beam is thus 
analogous to that between superbeam and reactor experiments m 

3 . 1 . Detection of the appearance signal 

The signal for the golden transition is a charged-current event (CC) with a muon in the 
final state. Reference I 2 n] studied the performance of a 440 kton fiducial water Cherenkov 
detector similar to Hyper-Kamiokande or the proposed by the UNO experiment |3()| . This 
analysis can be extended to different 7 ’s, using the same neutrino physics generator, detector 
simulation and reconstruction algorithms as described in uni, with realistic e-/x separation 
by pattern recognition, and the requirement of a delayed coincidence from muon decay. 

Figure [2 shows the reconstructed energy spectra of signal and background at maximum 
CERN-SPS 7 , for two different values of ^ 13 . Backgrounds are smaller for ®He than ^®Ne, 
and both neutrino and anti-neutrino backgrounds tend to cluster at low energies. Most of 
the background reconstructs below 500 MeV. 

The neutrino energy resolution depends strongly on the proportion of quasi-elastic (QE) 
and non quasi-elastic (non-QE) interactions in the signal. Neutrino energy is reconstructed 
assuming two-body, quasi-elastic kinematics, so contamination from non-QE events intro¬ 
duces a bias between the true and reconstructed energies. Figure |21 shows the fraction of QE 
and non-QE events passing the selection criteria. As expected the non-QE contamination is 
smaller for anti-neutrinos since the average beam energy is also smaller for the chosen 7 ’s. 

To properly include both detector resolution and non-QE contamination effects, a matrix 
describing the migration between true and reconstructed neutrino energies is constructed. 
Migration matrices are also computed for the backgrounds. Given the irreducible Fermi 
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Figure 2: Reconstructed energy for signal with 613 = 8 ° (solid) and 6*13 = 3° (dashed)) and 
background (dotted) at the maximum acceleration o/®He (left) and (right) ions at the 
CERN-SPS. The absolute normalization corresponds to one year. 




Figure 3; Quasi-elastie and non-quasielastic components in the p. appearanee signal for unit 
oscillation probability (the absolute normalization is arbitrary) at maximum CERN-SPS ac- 
eeleration o/®He (left) and ^®Ne (right). 

motion and muon threshold, the first energy bin extends from 0-500 MeV and bins of 250 MeV 
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Table 2: Migration matrices for appearance signal ) and backgrounds ) at different 
values of 'j. Each row and column of the matrices corresponds to a neutrino energy bin, as 
described in the text. 


width are used above it. For the high -7 Setup III, the first bin is discarded. 

Table 121 shows these migration matrices for 7 = 120 and 150 for ®He and ^®Ne. Three 
bins; 0-500, 500-750 and 750-1- MeV are used. The efficiencies are quite high (~ 30 — 50%) 
even when the background fraction is held below 10 “^. 

3.2. Detection of the disappearance signal 

For i^e —> t'e(7'e ^e) transitions, the signal is a CC interaction with an electron (positron) 

in the final state. In | 2 ^ this channel was included with a conservatively estimated 50% flat 
efficiency and negligible background. Since the energy resolution is also strongly affected by 
the non-QE contamination for this sample, this analysis is now refined to include the effect 
of migrations. While the background level for this large signal can be safely neglected in 
comparison to other systematic errors to be discussed later, a matrix of efficiencies should be 
used to account for the signal migrations. Table 01 shows these matrices for ®He and ^®Ne at 
various 7 ’s. Efficiencies are quite high, especially at lower energies where they reach 80-90%. 

3.3. Atmospheric background 

An important background for any accelerator-based experiment to control arises from atmo¬ 
spheric neutrinos. A detector like Super-Kamiokande will expect approximately 120 
interactions per kiloton-year (including the disappearance of into Vr). Of these, 32 atmo¬ 
spheric per kiloton-year pass all the selection cuts (one non-showering ring, accom- 

®Matrices with appropriate binning for other choices of 7 can be obtained from the anthors on request. 
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Table 3: Fractional migration matrices of the CC disappearance signal for different 

values ofj. 


7 

Selection 

Ejnax cut 

^min cut 

cos 61 cut 

120 

32 

19 

19 

15 

150 

32 

24 

24 

15 

350 

32 

30 

19 

5 


Table 4: Surviving atmospheric background per kton-year after cuts: on the high-energy 
end-point of the (5-beam neutrino spectrum (Emax), the low-energy tail (Emin) for setup III, 
and the lepton scattering angle (cos 9i), as described in the text. 


panied by a delayed coincidence from muon decay). The reconstructed spectrum of those 
events scaled by a factor 1/500 is shown in Figure 0 (solid line) alongside the signal for the 
three example setups to be considered later, namely, 7 = 120 (L = 130 km, dashed), 7 = 150 
{L = 300 km, dotted) and 7 = 350 L = 730 km, dashed-dotted), assuming ^13 = 1°. 

There are two additional handles to further reduce the atmospheric background. First, at 
a given 7 , we know the end-point of the signal spectrum, and there is no efficiency penalty 
for excluding events above the maximum beam energy. This cut obviously works best for 
lower -7 scenarios. Tabled shows the effect of the end-point cut for different 7 ’s. For higher 
7 , it is also helpful to set a lower energy cut. Requiring E > 500 MeV, for instance, is free 
for the highest 7 = 350 option, since this bin is not considered in the analysis anyway. 

Second, a directional cut is also possible, since the beam arrives from a specific, known 
direction but the atmospheric background is roughly isotropic. While the neutrino direction 
cannot be measured directly, it is increasingly correlated with the observable lepton direction 
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at high energies. Figure [5] illustrates this correlation for the three reference set-ups. Thus, a 
directional cut is more effective as 7 increases, but is never perfectly efficient. To compare 
the power of this cut for the different setups, we define it to achieve a 90% efficiency in all 
cases: cos 9i > 0.45 for 7 = 350, cos 9i > —0.3 for 7 = 150 and cos 9i > —0.5 for 7 = 120. 
The remaining atmospheric background for each setup is summarized in Table 01 Thanks to 
the directional cut, background rejection for the highest 7 is a factor three better than the 
alternative scenarios. 

Even with energy and directional cuts, 5 to 15 atmospheric background events per 
kiloton-year remain, compared to the expected intrinsic beam-induced detector background 
(mostly due to NC single-pion production) of 0(10“^) events. To reduce atmospheric con¬ 
tamination to a negligible level (say ten times below the intrinsic background) would require 
a rejection factor O(IO^) , although since the atmospheric background can be well measured 
a rejection factor 5-10 times less stringent is probably tolerable. 

This rejection factor can be achieved by timing of the parent ion bunches. It is esti¬ 
mated m that a rejection factor of 2 x 10^ is feasible with bunches 10 ns in length. Based 
on the present results, a less demanding scheme for the number of bunches and bunch length 
could be workable. 



Figure 4: Solid line: energy spectrum of atmospheric background per kiloton-year, 

scaled down by a factor 1/500. Dashed, dotted and dash-dotted lines: energy spectrum of 
signal events per kiloton-year for 7 = 120,150 and 350 assuming Ois = 1°. 
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3.4. Systematic errors 

Although a detailed analysis of all possible systematic errors is beyond the scope of this paper, 
we have included the two that will likely dominate. First, the uncertainty in the fiducial mass 
of the near and far detectors, which we estimate as a ±5% effect on the expected far-detector 
rate. Second, the uncertainty on the ratio of anti-neutrino/neutrino cross sections, which we 
assume a near detector can measure with an accuracy of ± 1 %^. 

To include these errors, two new parameters are added to the fits: A, the global normal¬ 
ization, and X, the relative normalization of anti-neutrino to neutrino rates. More precisely, 
if is the number of measured muon and electron events in the energy bin i for the anti¬ 
neutrino (-I-) or neutrino (—) beam, and ^"^’^(013, (5) is the expected number for some values 
of the unknown parameters (^ 13 , 5), then we minimize the following function: 


X^( 6 'i 3 ,( 5 , A,x) = 2 ^ <1 log 


n, 


AxNf 


+ ANj “ 


n 


I,— 


ANT 


iA-lf ix-lf 

“T O “T 


(y, 


erf 


(3.1) 


’^The calculable neutrino and anti-neutrino energy spectra of the /?-beam will facilitate cross-section mea¬ 
surements, compared to a traditional neutrino beam. 



cosei 


Figure 5: Cosine of the reeonstrueted neutrino-lepton seattering angle for three setups: 7 = 
120 (top), 7 = 150 (middle) and 7 = 350 (bottom). 
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where aA = 0.05 and ax = 0.01. The minimization in the parameters A and x for fixed values 
of 013 and 6 can be done analytically to leading order in the deviations A — 1 and x — 1, that 
is solving the linearized system: 


dA 


= 0 



(3.2) 


In what follows, sensitivity to the parameters (0i3,(5) will be quantified using 99% confi¬ 
dence regions for two degrees of freedom; that is, the curves satisfying: 


X iPlZt Amini Xmin) — 9.21. 


(3.3) 


4. Optimization of the CERN-SPS /3-beam 

The following sensitivity plots are used to optimize the physics performance of different (5- 
b earns: 

• Sensitivity to CP violation: region on the plane ( 0 i 3 ,h) where the phase 5 can be 
distinguished from both h = 0° and 5 = 180° for any best fit value of 0 i 3 , at 99% 
confidence level or better. 

• Sensitivity to ^ 13 : region on the plane ( 0 i 3 , 5) where the angle 0 i 3 can be distinguished 
from 013 = 0 for any best fit value of 5, at 99% confidence level or better. 

Unless otherwise specified, the following solar- and atmospheric-neutrino oscillation pa¬ 
rameters are assumed: 

Am ?2 = 8.2 X 10"%U^ 0 i 2 = 32° = 2.2 x IQ-^eU^ 023 = 45° (4.1) 

4.1. Optimal 7 for the CERN-Prejus baseline 

One frequently considered standard setup adopts the CERN-Frejus baseline L = 130 km 
and 7 = 60/100 for ®He/^®Ne [T^EI]. This setup appears to be far from optimal even if 
the baseline is kept fixed. As noted in [ISl, a higher -7 beam increases the event rate and 
allows the energy dependence of the signal to be analyzed. Taking the identical 7 for ®He and 
^®Ne, Figure ini shows the 7 -dependence of the 99% CL 6 and 0 i 3 sensitivity, as defined above. 
The stars indicate the values of the previously considered setup in nniEi], corresponding 
to 7 = 60/100. Clearly the CP-violation sensitivity is significantly better for larger 7 . For 
7 > 100 the sensitivity to CP violation and 0 i 3 changes rather slowly. This is not surprising, 
since increasing 7 at fixed baseline does not reduce the flux significantly at low energies (see 
Figure O, just as for a Neutrino Factory. In the absence of backgrounds, there is no penalty 
associated with higher 7 , although in practice, the non-negligible backgrounds result in a 
small decrease in 0 i 3 sensitivity at higher 7 , for some values of 6. 
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L=130 km 



'^ion 


Figure 6 : 'y-dependence of 99% confidence level 6 sensitivity at 613 = 8 ° (top) and 613 sen¬ 
sitivity (bottom) for 6 = +90° (solid) and 5 = —90° (dashed), assuming L = 130 km and 
76 He = 7 i 8 Ne- ^hc stars indicate the values for the 7 = 60/100 option in 11.51 Ufif . 


X 10'^ 



+(GeV) 

Figure 7: Energy spectra of (dashed) and Ve(solid) at L = 130 km for 7 = 100,120,150. 
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Figure 8 : Number of CC appearance candidates (from for unit oscillation probability, 

as a function of 7 , holding 7 /-L fixed. 


Although there is no unique optimal 7 within the wide range 7 = 100 — 150 when the 
baseline is fixed to L = 130 km, consider for illustration an intermediate 7 = 120 to define 
Setup I; a different choice of 7 > 100 will not make a significant difference. 

There appears to be no advantage to the asymmetric choice 7i8Ne/76He = 1-67. The 
asymmetric option is always comparable in sensitivity to a symmetric one with the smaller 
7 of the two, so a symmetric 7 conhguration is adopted for setup I. 

4.2. Optimal L for maximum ion acceleration 7 = 150 

As argued in [221, physics performance should improve with increasing 7 , if the baseline 
is correspondingly scaled to remain close to the atmospheric oscillation maximum, due to 
the (at least) linear increase in rate with 7 . This growth in sensitivity eventually saturates 
for a water detector, which becomes inefficient in reconstructing neutrino energies in the 
inelastic regime. Figure [HI where the number of CC appearance candidates selected (for unit 
oscillation probability) is plotted as a function of 7 (for 7 /L fixed), confirms this expectation. 
Saturation occurs for 7 ~ 400, above the maximum acceleration possible at the CERN-SPS, 
since the flux is still large in the quasi-elastic region (see Figure |7|). 

Fixing 7 to the CERN-SPS we next study the optimal baseline and how the symmetric 7 
setup compares with the asymmetric one. 

Figure ini shows tho |(5| ciiici ^13 SGiisitivitiGS cLS cL fuiictioii of the bsiSoliriG for — 150/150 
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99% CL sensitivity to CP violation, for 4.4 Mton y = 8 deg) 99% CL sensitivity to 813, for 4.4 Mton y 




Figure 9: Left: minimum value of\ 6 \ distinguishable from 0 and 180° at 99% CL (for 613 = 8 °) 
vs. baseline for j = 150/150 (red) and 7 = 150/250 (green). Right: minimum value of 613 
distinguishable from 0 at 99% CL (for 6 = 90° and —90° as shown). 


and the asymmetric case 7 = 150/250. The best CP sensitivity is achieved around L ~ 
300(350) km for symmetric(asymmetric) beams The baseline dependence of ^13 sensitivity 
leads to similar conclusions, although the importance of choosing the optimum baseline is 
more pronounced. A significant loss of 613 sensitivity results if the baseline is too short, as 
in Setup I. 

Setup II will hence be defined as 7 = 150/150 for L = 300 km. Similar results are expected 
for the asymmetric option 7 = 150/250 with slightly longer baseline. 


5. Comparison of the three setups 

From the results of the previous section, the default setups be compared are: 

• Setup I: 76 He = 7 i 8 Ne = 120 at L = 130 km 

• Setup II: 76 He = TisNe = 150 at L = 300 km 

• Setup III: 76 He = TisNe = 350 at L = 730 km 

For the highest 7 option, we have also checked that the symmetric and asymmetric options 
give comparable results. 

®The optimal baseline will obviously shift if is varied from the present best fit value: ±50 km for a 

change of one a. 
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99% CL sensitivity to CP vioiation, for 4.4 Mton y 99% CL sensitivity to 6 , 3 , for 4.4 Mton y 




Figure 10: Left: CP-violation exclusion plot at 99% CL for the three reference Setups I 
(dashed), II (dotted) and III (dashed-dotted) compared with the standard (solid) one of \15\ 
mif . Right: exclusion plot for ^13 at 99% CL with the same setups. The solar and atmospheric 
parameters are fixed to their present best fit values and the discrete ambiguities are assumed 
to be resolved. 


5.1. Intrinsic sensitivity to ^13 and 5 

Figure 11171 compares the CP-violation and @13 exclusion plots for the three setups assuming, 
for the moment, that the discrete ambiguities in sign(Am 23 ) and sign(cos 023 ) can be ignored 
because correct assignments have been made. Also included for reference is the previously 
considered setup from m- Although the highest 7 option of |25| remains best, the perfor¬ 
mance of Setup II is comparable. Even the sensitivity of the much-improved CERN-Frejus 
scenario in Setup I is considerable. Although only the range (—90°, 90°) is shown, to make 
it easier to read the y-scale, the region around 180° has a similar pattern. 

As explained in unj, differences between the setups arise due to sample size (which increases 
at least linearly with 7 ) and more robust energy reconstruction at higher energies (as Fermi 
motion becomes less important). 

Figure im shows typical fits for the three setups at several true values of 613 and 5. While 
both Setups II and III manage to resolve the intrinsic degeneracy essentially everywhere in 
the sensitivity range, this is not the case for Setup I; there (when the fake solution gets closer 
to and merges with the true one) the errors in 813 and 6 are sometimes strongly enhanced by 
the intrinsic degeneracy. This effect is not necessarily noticeable in the exclusion plot for CP 
violation. 
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Figure 11: Determination of { 613 , 6 ) at 99% CL for Setups III (thicker line), II (intermediate) 
and I (thinner line) and six different true values of the parameters indicated by the stars, 
assuming the correct sign(Am 23 ) and sign(cos 023 )- 


5.2. Effect of the eight-fold degeneracies 

By the time any /3-beam begins, it is probable that a number of uncertainties in the oscillation 
parameters besides ^13 and 6 will remain, in particular the discrete ambiguity in sign(Am 23 ) 
or the octant of 023- Both questions are theoretically important and the possibility of an¬ 
swering them with a /3-beam is attractive. These ambiguities are problematic, if they can’t 
be resolved, because they can bias the determination of the parameters ( 0 i 3 ,( 5 ), that is, the 
solutions surviving with the wrong assignment of the sign and/or the octant lie at different 
values of ^13 and d than the true ones. 

Generically, an eight-fold degeneracy of solutions appears when only the golden channel is 
measured and no energy dependence is available. There are two solutions in the absence of 
the discrete ambiguities, the true and the intrinsic one ini Each gets an false image for the 
wrong assignment of the sign m, for the octant da ED] and for both. 

As explained in [ 2 .‘ 1 | . the intrinsic solution and its three images are strongly dependent on 
the neutrino energy and therefore can be excluded, in principle, when the energy dependence 
of the oscillation signal is significant. On the other hand images of the true solution are 
energy independent and impossible to resolve unless there are additional measurements (e.g. 
disappearance measurements or the silver channel), or when there are significant matter 
effects. 
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Figure ca shows fits including the discrete ambiguities on the plane ( 0 i 3 , 5 ) for the three 
setups and different choices of the true ^13 and 6. In Setup I we generically find the full 
eight-fold degeneracy, while in Setups II and III the intrinsic solution and its images are 
typically excluded, thanks to the stronger energy dependence. 

Some general observations concerning these results include: 

• Presence of the intrinsic degenerate solution or its images as in Setup I is problematic, 
because it implies a signihcant increase in the measurement errors of ^13 and 6 (as 
shown in Figure ITT]) for some values of d). 

• When only the images of the true solution survive, as in Setups II and III, they interfere 
with the measurement of 613 and 5 by mapping the true solution to another region of 
parameter space. In vacuum [IH1I13]: 

Wrong-sign: 613 613 , d —> tt — 5 

Wrong-octant: 613 —> tan 023^*13 + 0 {Aml 2 ), sin5 ^ cot 623 sind 

Since these different regions occur for different choices of the discrete ambiguities they 
cannot overlap and one ends with a set of distinct measurements of 813 , 6 with different 
central values but similar errors (see the middle and right plots of Figure ini). 

• In vacuum, CP violating solutions are mapped into CP violating solutions, therefore 
the effects of degeneracies on the exclusion plot for CP violation are often small, even 
when degeneracies are a problem. In matter, on the other hand, 5 shifts in the fake 
solutions are enhanced by matter effects and for some central values of (^ 13 , <5) the fake 
solutions may move closer to the CP conserving lines than the true solution, resulting in 
an apparent loss of sensitivity to CP violation. This effect is visible in Figure IT^ where 
the fake-sign solution, which in vacuum should be located at ~ —140° for 5 = —40°, 
gets shifted towards the CP conserving line—180° for longer baselines where matter 
effects are larger. 

Figures El and m show the range of (^ 13 , 5), where the sign(Am 23 ) and sign(cos023) can 
be measured respectively. Asymmetric 7 options are also included, since there are some 
differences. As expected, sensitivity to the discrete ambiguities is better for large 613 and 
larger 7 . In Setup I there is essentially no sensitivity anywhere on the plane. 

Sensitivity to the discrete ambiguities and their bias in the determination of the parame¬ 
ters 013 and 6 could be significantly improved if data for any of the setups is combined with 
disappearance measurements, for instance in a superbeam experiment. This com¬ 
bination was recently studied in ESI for the standard /3-beam with significant improvement 
in sensitivity to the mass hierarchy, even without matter effects. A similar study for the 
setups considered here will be very interesting. One of the most important limitations of 
the /3-beam, compared to the superbeam or the Neutrino Factory, is its inability to measure 
the atmospheric parameters {92‘i, Arn^^) with precision. At the very least, information from 
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T2K phase-I should be included, since otherwise the uncertainty on these parameters will 
seriously compromise sensitivity to @13 and 6. Synergies in resolving degeneracies, between 
the /3-beam and T2K, should also be exploited. 

Another interesting observation is that atmospheric neutrinos can be measured in the same 
megaton detector considered here. A recent study m combining atmospheric data with T2K 
phase-II has found a large improvement in sensitivity of the latter to both discrete ambiguities 
when 013 is not too small (> 4°). This is surely an analysis that should be done and will be 
reported elsewhere. 


6 . Conclusions 

This paper has explored the physics potential of a CERN-SPS /3-beam, where ions can be 
accelerated to yeHe < 150 and yis^e < 250. The design of a /3-beam reaching this maximum 
7 is technically equivalent to the lower -7 option previously considered, for which a feasibility 
study already exists m A major improvement in sensitivity to 0 i 3 and 5 is achieved by 
increasing 7 . Even when the baseline is fixed to that of CERN-Frejus, sensitivity improves 
considerably if 7 > 100 and changes slowly as the 7 increases further to the limit of the SPS. 

An even more dramatic improvement is possible if the baseline is increased proportionally, 
so the first atmospheric oscillation maximum corresponds to the average neutrino energy, 
which occurs at L ~ 300 km. For large values of 0 i 3 this option is comparable in CP violation 
sensitivity to the optimal one in |25j at even higher 7 ~ 0(400), which would require a more 
powerful accelerator, such as the Tevatron or a refurbished SPS. In contrast, for small values 
of 013 the latter option is still significantly better. 

The main differences can be traced to increased event rate and the more significant energy 
dependence, which allows higher -7 options to resolve the intrinsic degeneracy. 

For discrete ambiguities, higher -7 also provides a window on the neutrino mass hierarchy 
and the octant of 023 , if non-maximal, relying on significant matter effects; the highest -7 
setup with 730 km baseline is therefore the only one with a significant sensitivity. 

In summary, if the existing CERN-SPS is the ion accelerator and the CERN-Frejus baseline 
is fixed, 7 should still be increased to a value greater than 100 , higher than considered 
in da ED- If an alternative site hosts a large underground laboratory near CERN, it will be 
profitable to exploit longer baselines L = 300 km. In any case, R&D effort to design /3-beams 
beyond the limit of the CERN-SPS appears justified, given the significant improvements in 
physics sensitivity they would allow. 
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Figure 12: Solutions for {6is,5) for the true‘ 2 ^lues: 6 = ±40° and Ois = 4° in (a) Setup I, 
(b) Setup II and (c) Setup III without discrete ambiguities, with the sign, octant and mixed 
ambiguities ordered from thicker to thinner-line contours. 





99% CL sensitivity to sign(Am^ 23 ) > O' V 09% CL sensitivity to signlAm^j) < 0, for 4.4 Mton y 




Figure 13: Region on the plane (^ 13 , 5) in which sign(Am 23 ) can he measured at 99% CL 
for 6»23 = 40.7° and positive (left) and negative (right) Symmetric and asymmetric 

beam options are shown for Setup II (300 km, solid and dashed, respectively) and Setup III 
(730 km, dotted and dash-dot). There is no sensitivity for Setup I. 


99% CL determination of 623 octant for Gjs = 40.7 deg (4.4 Mton y) 99% CL determination of 623 octant for 623 = 49.3 deg (4.4 Mton y) 




Figure 14: Region on the plane (0i3,(5) in which sign(cos023) can be measured at 99% CL 
for 023 = 40.7° (left) and 023 = 49.3° (right). Setups II (300 km, solid: symmetric, dashed: 
asymmetric) and III (730 km, dotted: symmetric, dash-dot: asymmetric) are shown. There 
is no sensitivity for Setup 1. 
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